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Abstract. We prove that the classical normal distribution is infinitely divisible with 
respect to the free additive convolution. We study the Voiculescu transform first by giving 
a survey of its combinatorial implications and then analytically, including a proof of free 
infinite divisibility. In fact we prove that a subfamily Askey-Wimp-Kerov distributions 
are freely infinitely divisible, of which the normal distribution is a special case. At the 
time of this writing this is only the third example known to us of a nontrivial distribution 
that is infinitely divisible with respect to both classical and free convolution, the others 
being the Cauchy distribution and the free 1/2-stable distribution. 



1. Introduction 

We will prove that the classical normal distribution is infinitely divisible with respect to 
free additive convolution. 

This fact might come as a surprise, since the classical Gaussian distribution has no special 
role in free probability theory. The first known explicit mentioning of that possibility to 
one of us was by Perez-Abreu at a meeting in Guanajuato in 2007. This conjecture had 
arisen out of joint work with Arizmendi [3]. 

Later when the last three of the present authors met in Bielefeld in the fall of 2008 they 
were led to reconsider this question in the context of investigations about general Brownian 
motions. We want to give in the following some kind of context for this. 

In [H] two of the present authors were introducing the class of generalized Brownian 
Motions (GBM), i.e., families of self-adjoint operators G{f) (/ G "H, for some real Hilbert 
space "H) and a state ip on the algebra generated by the G{f), given by 

7reP2(2n) (ij)e7r 

Here V2{'2n) denotes the set of pairings of 2n elements, and t is a weight function for such 
pairings. The concrete form of t determines the specific Brownian motion. 
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The most natural example for such a GBM is classical Brownian motion, where t(7r) = 1 
for all pairings vr; in this case one gets the normal law for G{f). The g-Brownian motion 
fits into this frame by putting tg{7i) = q^^^'^^ (where cr(7r) denotes the number of crossings 
of the pairing vr); in this case the law of the random variable G{f) is related with the theta 
function of Jacobi, and called g-Gaussian distribution 7^, see [21 [131 [12] • If the parameter 
q changes from -1 to 1, one gets an interpolation between the fermionic Brownian Motion 
{q = —1), the free Brownian Motion (g = 0), and the classical Brownian Motion (g = 1). 

In [11], the model of the free product of classical Brownian motions resulted in a new 
class of GBM, with the function t given by 

(1.1) Un) = s^^W, 

where cc{7c) is the number of connected components of the pairing vr. Here s has to be 
bigger than 1. This contains as a special case the result: The 2?7.-th moment of the free 
additive power of the normal law 71 is given as follows: 

7reP2(2n) 

for s > 1 . 

In the light of earlier examples where similar combinatorial identities could be extended 
beyond their primary domain of applicability (see, e.g., [10]), it was natural to ask whether 
this relation could also make any sense for s < 1. So, in this context, a natural problem 
is whether the sequence on the right side of the above formula is a moment sequence 
for all s > 0? This question is equivalent to the free infinite divisibility of the normal 
law! One can check easily that the corresponding more general question on generalized 
Brownian motions (i.e., whether the tg from equation f 1 1.11) is still positive for s < 1) has 
a negative answer. From this point of view, the free infinite divisibility of the classical 
Gauss seemed quite unlikely. However, numerical evidence suggested the validity of that 
conjecture. In this paper we will give an analytical proof for this conjecture. We want 
to point out that it still remains somehow a mystery whether the ffl-infinite divisibility of 
the Gauss distribution is a singular result or whether there is a more conceptual broader 
theory behind this. 

Another example of this phenomenon was found in ^35j, namely that the 1/2- free stable 
law [7\ is also classically infinitely divisible, being a /3- distribution of the second kind with 
density (4x — l)^/^/a;^. 

1.1. Related Questions. One way to describe certain probability distributions on M is 
by specifying their orthogonal polynomials. The orthogonal polynomials of the classical 
Gaussian distribution are the so-called Hermite polynomials [27j. In |4J, Askey and Wimp 
describe a family of deformations, indexed by c G (—1, +c>o), of the Hermite polynomials, 
called the associated Hermite polynomials. These polynomials are orthogonal with respect 
to a family of probability measures {/ic: c G (—1, +00)}, which can be described in terms 
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of a continued fraction expansion of their Cauchy-Stieltjes transform as 




1 



c + 1 



z — 



c + 2 



z — 



z — 



The Cauchy-Stieltjes transform of a measure \i on the real line is defined by 



z G C+; 



when /i is a positive measure, this function maps C"*" into the lower half-plane. For details 
we refer to the excellent book of Akhieser [1]. For c = 0, we have /iq = 7, the normal 
distribution |27j, and one can easily check that one can extend this family continuously to 
c = — 1 by letting = ^o, the probability giving mass one to {0}. This family, introduced 
in [1], plays an important role in [27]; we shall call its members the Askey-Wimp-Kerov 
distributions. It will turn out from our proof that {iJ,c'- c G [—1,0]} are freely infinitely 
divisible. Numerical computations show that for several values of c > 0, /ic is not freely 
infinitely divisible. Numerical evidence seems also to indicate that /^c is classicaly infinitely 
divisible only when c = Oorc = — 1. 

An interesting interpolation between the normal and the semicircle law was constructed 
by Bryc, Dembo and Jiang [IB] and further investigated by Buchholz [19]. This leads to a 
generalized Brownian motion, given by a weight function (0 < 6 < 1) 



where for a pairing vr of 2n elements we denote by /i(vr) the number of connected components 
which have only one block with 2 elements. It was calculated in [TSl [12] , that the measures 
satisfy 



where Db is the dilation of a measure by parameter 6 > 0. With the results of the 
present paper it is immediate that these measures are infinitely divisible with respect to 
free convolution. However a short calculation shows that these measures are not infinitely 
divisible with respect to classical convolution unless 6 = 1. 

One other tempting example is to consider the distribution of N x N Gaussian random 
matrices. For = 1, this is the classical Gauss distribution, whereas for A^ — t- oo it 
converges to the semicircle distribution. Both of them are infinitely divisible in the free 
sense (the prior by our main result here, and the latter because the semicircle is the 
limit in the free central limit theorem). So one might conjecture that the interpolating 
distributions, for integer 1 < A^ < oo, are also freely infinitely divisible. However, numerical 
calculations of the first few Jacobi coefficients of the corresponding moments, using the 
Harer-Zagier recurrence, show readily that this is not the case. 



Z/;,2 = Db'Jl ffl I^l_b70, 
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One may also ask the "opposite" question, whether the Wigner distribution is infinitely 
divisible with respect to classical convolution. However this is impossible because any non- 
trivial classically infinitely divisible measure has unbounded support, see ^21 Proposition 
2.3]. 

For the same reason the distributions whose density is a power of the Wigner density 
are not classically infinitely divisible. It is however an open question whether the latter are 
freely infinitely divisible. Numerical evidence points to a positive answer to this question. 
This would provide another proof that the normal law is freely infinitely divisible. See [3] 
for a survey on these questions. 

In the next section we will consider the combinatorial aspects of the free infinite divisibil- 
ity of the classical Gaussian distribution; in particular, we will give some new combinatorial 
interpretations for the free cumulants of the Gaussian distribution. In Section 3, we will 
then give an analytical proof of our free infinite divisibility result. 

Acknowledgements The last three authors would like to thank Prof. Gotze for the kind 
invitations and hospitality at SFB 701, Bielefeld. STB would also like to thank Michael 
Anshelevich for many useful discussions, especially regarding the Riccati equation. 

2. Combinatorial considerations 

2.1. Partitions. First we review a few properties of set partitions which will be needed 
below. As usual, set partitions will be depicted by diagrams like the ones shown in Fig- 
ure [2Tl 

Definition 2.1. A partition of [n] = {1,2, ... ,n} is called 

(1) connected if no proper subinterval of [n] is a union of blocks; this means that any 
diagram depicting the partition is a connected graph. 

(2) irreducible if 1 and n are in the same connected component, i.e., there is only one 
outer block. 

(3) noncrossing if its blocks do not intersect in their graphical representation, i.e., if 
there are no two distinct blocks Bi and B2 and elements a, c G -Bi and h,d & B2 
s.t. a < b < c < d. Equivalently one could say that a partition is noncrossing if 
each of its connected components consists of exactly one block. 

Typical examples of these types of partitions are shown in Fig. 12. 1[ 

We denote the lattice of partitions of [n] by 7^, the irreducible partitions by T^^^ and 
the order ideal of connected partitions by V^°^^; the lattice of noncrossing partitions will 
be denoted by iVC„, and the sublattice of irreducible noncrossing partitions by NC^^. 
Finally, let us denote by X„ the lattice of interval partitions, i.e. the lattice of partitions 
consisting entirely of intervals. 



connected irreducible noncrossmg 



Figure 1. Typical partitions 
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2.2. Cumulants. Cumulants linearize convolution of probability measures coming from 
various notions of independence. 

Definition 2.2. A non- commutative probability space is a pair {A, 0) of a (complex) unital 
algebra A and a unital linear functional 0. The elements of A are called (non-commutative) 
random variables. 

Given a notion of independence, convolution is defined as follows. Let a and b be 
"independent" random variables, then the convolution of the distributions of a and b is 
defined to be the distribution of the sum a+b. In all the examples below, the distribution of 
the sum of "independent" random variables only depends on the individual distributions of 
the summands and therefore convolution is well defined on the level of probability measures. 
Moreover, the n-th moment m„(a + 6) is a polynomial function of the moments of a and 
b of order less or equal to n. For our purposes it is sufficient to axiomatize cumulants as 
follows. 

Definition 2.3. Given a notion of independence on a noncommutative probability space 
{A,(f)), a sequence of maps a i->- kn{a), n — 1,2, . . . is called a cumulant sequence if it 
satisfies the following properties 

(1) kn{a) is a polynomial in the first n moments of a with leading term m„(a). This 
ensures that conversely the moments can be recovered from the cumulants. 

(2) homogeneity: A;„(Aa) = A"A;„(a). 

(3) additivity: if a and b are "independent" random variables, then kn{a + b) = kn{a) + 
kn{h). 

Mobius inversion on the lattice of partitions plays a crucial role in the combinatorial 
approach to cumulants. We need three kinds of cumulants here, corresponding to clas- 
sical, free and boolean independence, which involve the three lattices of set partitions, 
noncrossing partitions and interval partitions, respectively. Let X be a random variable 
with distribution and moments m„ = m„(X) = j x"' dip{x) 



2.3. Classical cumulants. Let 



/oo 
e^-'dij{x) = Y, 
n 



OO 

nl 

n=0 



be the formal Laplace transform (or exponential moment generating function). Taking the 
formal logarithm we can write this series as 



where 



oo 

K{z) = E ^ ^" 

n=l 



is the cumulant generating function and the numbers «;„ are called the (classical) cumulants 
of the random variable X. 
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If for a partition vr = {tti, 7r2, . . . , VTp} we put m^r = ?^|7ri|'^|7r2r ' ' "^^^1 and k^, = 
' ' ' i^\tt |; then we can express moments and cumulants mutually using the Mobius 
function /i on the partition lattice n„ as follows: 

a<.7V (T<.1T 

For example, the standard Gaussian distribution 7 = A^(0, 1) has cumulants 



1 n = 2 
n^2 



It follows that the even moments = of the standard gaussian distribution count 
the number of pairings of a set with the corresponding number of elements. 

2.4. Free Cumulants. Free cumulants were introduced by Speicher |38] in his combi- 
natorial approach to Voiculescu's free probability theory. Given our random variable X, 
let 

00 

(2.1) M{z) = 1 + J2^nz'' 

n=l 

be its ordinary moment generating function. Define a formal power series 

00 

n=l 

implicitly by the equation 

C{z) = C{zM{z)). 

Then the coefficients c„ are called the free or non- crossing cumulants. The latter name 
stems from the fact that combinatorially these cumulants are obtained by Mobius inversion 
on the lattice of non-crossing partitions: 

(2.2) = ^ = ^ m^/iArc(cr, vr) 

a£NCn aeNCn 
(T<7r (7<7r 

2.5. Boolean cumulants. Boolean cumulants linearize boolean convolution Let 
again M{z) be the ordinary moment generating function of a random variable X defined 
by (12. ip . It can be written as 

where 



n=l 
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and the coefficients are called boolean cumulants. Combinatorially the connection between 
moments and boolean cumulants is described by Mobius inversion on the lattice of interval 
partitions: 

(2.3) m^=^^K = ^ m^iix{o;n) 

(T<7r a<.n 

The connection between these kinds of cumulants is provided by the following theorem 
(see also |T1] for the case of pairings). 



Theorem 2.4 ([29j). Let {rrin) be a (formal) moment sequence with classical cumulants 
Kn- Then the free cumulants of mn are equal to 

(2.4) c„ = ^ 

the boolean cumulants are equal to 

2.6. The normal law and pair partitions. Let (i7(t) = '^^~~dt be the standard 
normal (classical Gaussian) distribution. Its moments are given by 



ruk := / t^d-f{t) 




= (A;-l)(A;-3)(A;-5)---5-3-l, even 

k odd 

A more combinatorial description of this is that the moments count all pairings, i.e., 

rrik = m{k). 

From Theorem 12.41 it follows then that the free cumulants Cn of 7 are given by the number 
of connected (or irreducible) pairings, 

Cn = #{7r G 7^(n) I 7r is connected}. 

The question whether 7 is infinitely divisible in the free sense is equivalent to the question 
whether the sequence (c„)„gN is conditionally positive, which is the same as the question 
whether the shifted sequence (s„)„>o, where Sn '■= Cn+2, is positive definite, i.e., the moment 
sequence of some measure. See Section 13 of [33] for more details on this (note that there 
only compactly supported measures are considered, but the theory also extends to measures 
which have a uniquely solvable moment problem). 

The first few values of the free cumulants of the Gaussian distribution are 

02 = 1, C4 = 1, C6 = 4, 08 = 27, Clo = 248, 012 = 2830, . . . 

This sequence of the numbers of irreducible diagrams of 2n nodes has been well-studied 
from a combinatorial point of view, see, e.g., [HI HQ]; it appears for example as sequence 
A000699 in Sloane's Encyclopedia of Integer Sequences [37]. For a recent bibliography of 
this sequence see [28] where it is shown that the sequence is not holonomic, i.e., it does not 
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satisfy a linear recurrence with polynomial coefficients. However, positivity questions for 
this sequence have never been considered. It might be interesting to point out that these 
numbers appear also in the perturbation expansion in quantum field theory for the spinor 
case in 4 spacetime dimensions, see [15j and in renormalization of quantum electrodynam- 
ics, see [16]; however, due to the cryptic style of the mentioned papers the meaning of this 
remains quite mysterious for the present authors. 

2.7. A recursive formula. A main result on the numbers of irreducible diagrams is the 
following recursion formula due to Riordan 

n-l 

C2n = - 1) ^ C2iC2(„-i), 

a simple bijective proof of which can be found in [31]. In terms of the shifted sequence 
(ySn)n>o this reads 

n—l 

(2.5) S2„ = n ^ S2iS2{„-i-i). 

Note the similarity to the standard recursion of the Catalan numbers (just remove the 
factor n before the sum). An equivalent formulation is 

n-l 

(2.6) S2„ = ^(2i + l)s2jS2(n-i-l)- 

1=0 

Thus the question is whether the sequence {sk)k>o ^ defined by S2n+i = (n G N) and 
by either of the recursions (12. 5p or (12. 6p and Sg = 1 - is the moment sequence of some 
measure. Since Sq = 1, this measure must necessarily be a probability measure. The most 
direct way to prove this would be to find a self adjoint operator which has these numbers 
Sn as moments. 

There are some immediate combinatorial interpretations of the above recursions. For 
example, the recursion (12. 6p yields 

Here we are summing over all non-crossing pairings of 2n elements and the contribution of 
a pairing vr is given by a product over the blocks of vr, each block contributing the number 
ip(y) of its inner points plus one. These inner points have also been counted in [il| [23] in 
different contexts. 

2.8. Tree factorials. The recursion (12. 5p . on the other hand, can be interpreted in terms 
of planar rooted binary trees. These are planar rooted trees such that each vertex at most 
2 successors, called children. A vertex without successors is called a leaf. Denote the set 
of such trees with n vertices by VTZBTn- The number of these trees is the n-th Catalan 
number. The tree-factorial is defined as follows. For n = there is only one binary tree 
(the empty tree), whose factorial is defined to be 1. Let t be a binary tree with n > 
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vertices. Then t can be decomposed into its root vertex, a left branch ti with k vertices 
and a right branch ^2 with n — 1 — k vertices and we define 

t\ =n-ti\t2\ 

Then we have the following identity. 
Proposition 2.5. 

(2.7) S2n= Yl 

Indeed using the above decomposition it is easy to see that the numbers on the right 
hand side also satisfy the recursion (12. 5p . For more information on tree factorials see, e.g., 
PT| Section 2] and [301 Section 2] and section I2.1UI below. 

Note that these interpretations are canonical for the shifted sequence (s^) and not for 
the original sequence (c„); for example, Cg = 27 is the number of irreducible pairings of 8 
points, but = 27 is given in terms of non-crossing pairings of 6 points or, equivalently, 
in terms of planar binary trees with 3 = 6/2 nodes. 

2.9. Two Markov chains. 

2.9.1. MTR on binary search trees. The tree factorial appears in the stationary distribution 
of the move-to-root Markov chain on binary trees [22] . Binary trees are used in computer 
science to arrange data such that it can be accessed using binary search. To reduce search 
time, every time an entry is searched it is moved to the root of the tree by repeating the 
so called simple exchange shown in the following picture 




until the root position is reached. Choosing a vertex randomly (each with probability 1/n), 
this induces a Markov chain on the state space VTZBTn- By Perron- Frobenius theory there 
is a unique stationary distribution vr for this Markov chain and it is shown in [22] that it is 
given by Tx{t) = l/t\. Equivalently, it describes the distribution of a randomly grown tree. 

2.9.2. The Naimi-Trehel algorithm on planar rooted trees. The tree factorial also appears 
in the so-called Naimi-Trehel algorithm [43^ i32j. This is a queuing model based on yet 
another Catalan family, namely planar rooted trees VTZT . It solves a scheduling problem 
for n clients (e.g., computers) who access some resource (e.g., a printer) which can serve at 
most one client at a time. In order to reduce the number of messages needed to schedule 
the printer jobs, the queue is arranged as a planar rooted tree and each time a request is 
sent, the queue is rearranged. The average number of messages is then a certain statistic 
on these trees. This can be modeled as a Markov chain on labeled rooted trees where at 
each step a random client sends a request and the tree is transformed accordingly. In the 
end only the shape of the tree matters and it suffices to consider the corresponding Markov 
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chain on the unlabeled planar rooted trees. By means of bijection this can be transformed 
into a Markov chain on Dyck paths where we have the following algebraic rule for the 
transition probabilities [32] . 

Let us consider words in the two letter alphabet {x, x*} where x is an upstep or NE step 
and X* is a downstep or SE step. A Dyck word is a word in x and x* such that each left 
subword contains not more downsteps than upsteps and the whole word contains an equal 
number of up- and downsteps. Dyck words can be visualized by Dyck paths, see fig. [2l 
These are lattice paths which do not descend below the x-axis. We denote by 1 the Dyck 

.-v^ .-v-v-s 

Figure 2. Dyck paths of length 6. 

word of length and Dn the set of Dyck words of length 2n. 

The recursive structure of rooted planar binary trees has a counterpart in the unique 
decomposition of a Dyck word w as concatenation 

w = uxvx* 

such that both u and v are again Dyck words. Using this recursive structure we get a 
natural bijection a from planar rooted binary trees to Dyck words by setting recursively 
a{t) = a(ti)xa{t2)x* if t has left and right subtrees ti and t2- Under this bijection the tree 
factorial on Dyck words can be recursively computed as 

1! = 1 {uxvx*)\ = n ■ u\v\ 

Dyck words form a monoid with the concatenation product and following [32] we recursively 
define a linear operator on the monoid algebra of formal linear combinations of Dyck words 
by letting 

= w + ^{w) 

where 

z/(l) = u{uxvx*) = u{u) ® {xvx*) + fi{v)xux* 
and the operation ® is defined as 

uxvx* & w = uxvwx* 
For example, using these rules we have 

M/^) ^/^. 

Then it is easy to see by induction that n maps a Dyck path w of length 2n to a linear 
combination of Dyck paths of the same length and that the coefficients are nonnegative 
integers which sum up to n + 1. We interpret the matrix representation A of this map as a 
weighted adjacency matrix and obtain a digraph with vertex set Dn- Dividing the matrix 
by n + 1 we obtain a stochastic matrix P = j^A. It was shown in [32] that this is exactly 
the transition matrix of the Naimi-Trehel Markov chain discussed above. Figures [3] and H] 
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show the graphs for n = 2, 3 and 4. Again there is a unique stationary distribution tt for 




Figure 4. Naimi-Trehel graph of order 4. 



this Markov chain. It is shown in [32] that niw) = ^ where w\ is the tree factorial defined 

! i ^ ' ui! 

above. 
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2.9.3. A probabilistic interpretation. For both Markov chains discussed above our sequence 
appears as 



(2.8) S2n= V ^ = Ve^T^ = C„ET 



where by standard Markov chain theory 
(2.9) E^T, 



is the expected time of a random walker starting in w to come back to w for the first time 
and ET is the expected return time of a random walker starting at a randomly chosen 
state w, each chosen with probability 1/C„ (Catalan number). Although this setting is 
very close to the Lindstrom-Gessel-Viennot theory of determinants, we did not manage to 
exploit it for a combinatorial proof of our theorem. 

2.10. Loday-Ronco Hopf algebra. Hopf algebras of trees have enjoyed increasing inter- 
est recently in renormalization theory and noncommutative geometry [21] and pure algebra 
( "dendriform algebras" ) |30] . Hopf algebras of labeled trees have been studied by Foissy 
[22125]. 

We provide here yet another Hopf algebra on labeled trees whose Hilbert series is related 
to our problem. 

Definition 2.6. Let t be a planar rooted binary tree. A labeling of t is a function from the 
vertices of t to the integers. A labeling is called anti-increasing if the labels are distinct 
and for every vertex v of t, the labels of the left subtree (with root v) are strictly smaller 
than the labels of the right subtree. In other words, if we interpret the tree as the Hasse 
diagram of a poset, every antichain has increasing labels. Two trees with anti-increasing 
labelings are called equivalent if the induce the same linear order on the vertices. The 
equivalence classes are called anti-increasingly ordered trees. 

Proposition 2.7. Let t be a planar rooted binary tree with n vertices. Then the tree fac- 
torial t\ counts the number of anti-increasing orderings or equivalently the anti-increasing 
labelings with different numbers 1,2, ... ,n. 

The proof is a simple induction using the recursive definition of the tree factorial. 
Figure [5] shows an example of a anti-increasing tree. The formal linear combinations of 




Figure 5. An anti-increasingly ordered tree 



labeled planar rooted binary trees form a graded vector space, the grading being given by 
the number of vertices of the trees. It is then straightforward to generalize the coproduct 
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of Loday and Ronco [30] to labeled trees as follows |2ll |25]. Let s and t be labeled binary- 
trees. We define a new labeled binary tree s t by grafting them on a new root with label 
k: 

Contrary to the Butcher-Connes-Kreimer Hopf algebra of rooted trees, binary trees 
cannot be grafted from forests. Therefore the left part of the Loday-Ronco coproduct does 
not consist of forests, but rather a certain noncommutative "dendriform" product of binary 
trees. The product of two labeled binary trees s = si S2 and t = ti V; t2 is recursively 
defined as 

(2.10) S*t = SiVk{s2*t) + {s*ti)Vit2 

with the convention that for the empty tree | the product is 

\*t = t*\=t. 

This operation is associative and the coproduct of t = u f is recursively defined as 

A{t) = A{u) &k A{v)+t^ I 
where in Sweedler's notation we define for A{u) = J^'^W ^"^(2) and A{v) = ^"^(2') 

A(n) ®fc A{v) = ^ * ^(1')) ® («(2) Vfc V(2'))- 

Proposition 2.8. (1) Let s and t be anti-increasingly labeled trees such that the labels 
of s are smaller than the labels oft. Then s*t is a sum of anti-increasingly labeled 
trees. 

(2) Let t be a anti-increasingly labeled tree. Then all terms of A{t) contain anti- 
increasingly labeled trees only. 

Proof. (1) By induction, let s = Si Vfc S2 and t = ti V; t2, then s * t = Si \/ k {s2 * t) -\- 
{s *ti) V; ^2- By induction hypothesis, the monomials of both S2 *t and s * ti are 
anti-increasingly labeled, the labels of si are smaller than the labels of S2 * t and 
the labels of s * ti are smaller than the labels of ^2- 
(2) Assume that t = uVkV has anti-increasing labels and its children have coproducts 
^('") = Yl ^(1) ® ^(2) and A(f ) = f (2/), then by the preceding calculation 

each term * t>(i/) is a sum of anti-increasingly labeled trees and each U(2) V(^2') 
has anti-increasing labels as well, therefore 

A{t) = ^ ^{U(l) * f (1')) ® (M(2) Vfc f (2')) + t ® I 

only contains anti-increasingly labeled terms. 

□ 

It follows that the anti-increasingly labeled trees form a graded sub-coalgebra, but not a 
subalgebra because the product of anti-increasingly labeled trees does not consist of anti- 
increasingly labeled trees in general. Passing to anti-increasingly ordered trees we obtain 
a Hopf algebra. 
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Corollary 2.9. The anti-increasingly ordered trees span a Hopf algebra whose Hilbert series 
is the generating series of our sequence (s„). 

Proof. We have seen that anti-increasingly labeled trees form a coalgebra. Define the 
product of anti-increasingly ordered trees s and t as follows: Put arbitrary anti-increasing 
labelings on s and t such that the labels of s are smaller than the labels of t. Then 
s * t consists of anti-increasingly labeled trees and the corresponding anti-increasingly 
ordered equivalence classes do not depend on the choice of the labelings for s and t. Define 
the coproduct of an anti-increasingly ordered tree by choosing an ant i- increasing labeling 
compatible with the ant i- increasing order, compute the coproduct of the obtained anti- 
increasingly labeled tree and replace the anti-increasingly labeled terms of the result by 
the corresponding anti-increasingly ordered trees. Again the choice of the initial anti- 
increasing labeling has no infiuence on the final result; moreover, the obtained coalgebra is 
graded and connected, i.e., the first homogeneous component is one-dimensional, and the 
existence of the antipode follows by standard Hopf algebra theory. □ 

Some examples: 

A(|) = |®| 
A(Q)) = I ®® + ®® I 

A(^) = |®^ + ®®® + ^®| 
A(Y) = |®Y + ®«\ + ®®/ + /^® + %®® + Y^| 



2.11. Charge Hopf algebra. There is another coproduct ("charge Hopf algebra" of 
Brouder and Frabetti) in [17J, also mentioned in [26l sec. 2.6] which is more asymmet- 
ric than the Loday-Ronco Hopf algebra. First we define an associative multiplication s/t 
by putting \/s = s/\ = s and otherwise grafting s onto the leftmost leaf of t. Recursively, 
s/(ti Vt2) = (-5/^1) Vt2- This makes sense for labeled trees as well (just keeping the labels) 
and for LR-ordered trees we define s/t by shifting the labels of t such that all labels of s are 
less than the labels of t before carrying out the product. Then the algebra is generated by 
all elements of the form | t which we denote by Vk{t). The Brouder- Frabetti coproduct 
is defined recursively as 

A^(|) = |®| 

A^(®) = ®® I + I ®® 
A^Vkis t)) = V,{s Vit)(^\+td® Vk{A^s)/{A^Vi{t)) -Vi{t)(^\) 
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Again the monotone trees form a subalgebra. Examples: 




2.12. Conclusion. These interpretations of the sequence s„ invite to Fock space like con- 
structions with corresponding creation and annihilation operators, with the hope that the 
sum of creation and annihilation operator would have the s„ as its moments. However, 
we were not able to implement this idea successfully. The involved inner products usually 
lacked positivity. 

In the next section we will give an analytic proof of the positive definiteness of the 
sequence Sn- Essentially, it will consist in showing that the generating power series 

oo ^ 
n=0 

is actually the Cauchy-transform of a probability measure. We are not able to determine 
this measure directly, but we will show its existence. Note that the recursion (12. 6p implies, 
at least formally, for the equation 

This equation will be the starting point of our investigations in the next section and we 
will show that it allows an extension of (p as an analytic map from the upper to the lower 
complex half-plane, which is a characterizing property for Cauchy transforms. 

Since the Sn grow of the same order as the moments of the Gauss, the above formal series 
has no non-trivial radius of convergence. However, it does determine uniquely an analytic 
map on some truncated cone at infinity; our proof will show that this map extends to the 
upper half-plane. 



3. Analytic proof of the theorem 

In this section we shall give an analytic proof of the free infinite divisibility of the 
classical normal distribution. We shall obtain the free infinite divisibility of the classical 
Gaussian as a limiting case of a more general family of freely infinitely divisible distributions 
with noncompact support, namely the Askey-Wimp-Kerov distributions, defined in the 
introduction. A certain sub-family of the Askey-Wimp-Kerov distributions appears in 
[TT] . Recall that the distributions {fie'- c G (— l,-|-oo)} are determined by the continuous 
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fraction expansion of their Cauchy-Stieltjes transforms: 
(3.1) G,M = ^ 



z 



c + 2 
z 



z — ■ ■ 

In this section, we shall prove the following theorem: 

Theorem 3.1. With the notations from above, the probability measures Hc are freely in- 
finitely divisible for all c G [—1, 0]. 

To prove this result, we shall use the well-known characterization of free infinite divisi- 
bility provided by Bercovici and Voiculescu. Recall [8] that 

Theorem 3.2. A Borel probability measure fi on the real line is ^-infinitely divisible if 
and only if its Voiculescu transform (f)fj.{z) extends to an analytic function cj)^: C+ — t- C^. 

We remind the reader that the Voiculescu transform of a probability measure is defined 
by the equality (f)^{l/z) = Gj^'^{z) — 1/z, for z in some Stolz angle in the lower half-plane, 
with vertex at zero|3 For more details and important properties of this transform, we refer 
to [8]. In particular, this theorem guarantees that taking weak limits preserves free infinite 
divisibility. Our main source for the analysis of the function will be Kerov's work [27] 
and the paper |1] of Askey and Wimp. It is shown there that /ic is absolutely continuous 
with respect to the Lebesgue measure on the real line and 

3./) ; = , ; —, M G M, 

where for c > 0, 



,2 



DJz) = —— I e~'''-'^x'-Ux. 



Tic) 

This integral representation does not hold for c G (—1, 0], but the function D^diu) is still 
well defined for c G (—1,0), according to a theorem of Askey and Wimp [271 Theorem 
8.2.2]. Moreover, as remarked in [31 Section 4], the function D_c{z) is an entire function 
of both c and z. 

Remark 3.3. (1) We observe that, for fixed c > —1, equation (13. 2p together with the 
analyticity of Z) in c and z, guarantees that D^d^u) ^ for all m G R. Indeed, a 
zero of M 3 M D^d^u) G C would have to be of order at least one, hence the 
function |^ \iu)\'^ would not be integrable around that particular zero. 

(2) On the other hand, analyticity of D_c alone guarantees that the density '^^f^^ is 
strictly positive everywhere on the real line. 



""^Thc name "Voiculescu transform" for cj)^^ has been introduced in a later paper [7 of Bercovici and 
Pata, not in [S]. 
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In the same paper of Kerov [27j it is shown that G^^{z) satisfies the Riccati equation 

(3.3) G'^^iz)=cG,M'-zG,M + h ^eC+. 

It is also shown that this expression is equivalent, via the substitution 

which holds for any c 7^ 0, to the second order linear differential equation 

(3.4) ip"{z) + zip\z) + ap{z) =0, zE C+. 

(The function (p does depend on c.) The density ^ is analytic around zero, so, according 
to [HI Lemma 2.11], it follows that the function has an analytic extension to a small 
enough neighbourhood of zero. Clearly G^^(O) 7^ 0, so has an analytic extension around 
zero. Using fl3.4p . we obtain a convergent power series expansion for 99, namely f{z) = 
Y^'^=oCnZ'^, where Cn+2 = - („+i)"("+2) C«. We conclude that 

c{c + 2){c + A)---{c + 2k) 
C2k = (-1) CO, 

and 

,,, (c+l)(c + 3)---(c + 2fc-l) 
C2m - (-1) ^^^^y^ ci. 

Trivially, these coefficients provide a power series with an infinite radius of convergence. 
Thus, (p is an entire function. We immediately conclude that G^^ extends to a meromorphic 
function defined on all of C, whose poles coincide with the zeros of ip. We shall denote this 
extension also by G^^. (It may be worth noting that ip is entire for any c E C, and that 
when c — )■ —1, we have p{z) — )• z, hence /ic tends in the weak topology to the Dirac point 
mass at zero, 6q.) We shall denote 

(3.5) F,^{z) = -\-, zeC. 

Again, from the above it is clear that is meromorphic, its poles coinciding with the 
critical points of p. It satisfies the differential equation 

(3.6) F'^S^) = -F,^izy + zF,Sz) - z E C. 

(This follows easily when we divide by —G^^{zy in equation (13. 3p .) Dividing equation 
dMl) by Ff^^z) gives 

"^'^^^'^ - {z - F,Sz)) ' 



Proof, (of Theorem 13. 1[ ) We will prove our theorem by arguing that, essentially, the 
meromorphic function F^^ maps some simply connected domain containing the upper half- 
plane bijectively onto This, together with Theorem 13. 2[ will allow us to conclude. 
For further reference, we shall split the proof in a succession of remarks and lemmas. For 
the rest of the proof, we shall fix a c G (—1,0) (open interval!) Let us start by studying 
the critical points of F^^. It is known [1, Chapter III] that '^F^^z) > '^z for all z E C"*". 
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Moreover, if Ff^^{z) G while z C^, then it is clear that '<s{z — F^^{z)) < 0. In addition, 
as c < 0, (^—y^^ < whenever F^^{z) G C+. This, together with Remark [X^ implies 

Remark 3.4. If QF^^{z) > 0, then F'^^{z) ^ 0. In particular, F^^ has no critical points 
anywhere in the closure of the upper half-plane. 

It will be of use to also study the behaviour of F^^{iy), y G M. 

Remark 3.5. The continued fraction expansion of indicates that the probability 
measure /ic is symmetric with respect to the origin, so in particular G^^(i[0, +oo)) C 
i{—oo, 0]. This, of course, together with the meromorphicity of and G^^, requires that 
G^^{m.),F^^{m) CtRu{oo}. 

Part of the following lemma will not be used directly in the proof of the current theorem, 
but we find that it is nevertheless worth mentioning these results. 

Lemma 3.6. (1) There exists — oo < qo < so that F^^ maps i[go,+oo) bijectively 
onto i[0, +oo). 

(2) oo, go)) is a bounded subset o/z(— oo,0). In particular, F^^ has no poles on 
the imaginary line. 

(3) We have '^F^Xw) > V for all y eR. 

(4) The function F^^ has a unique (simple) critical point is in z(— oo, — 2-^/^). In 
addition, liiRy^^^ F'^^{iy) = \imy^_^i~^F^^{iy) = 0. 

(5) F'Jty)<lfor allyER. 

Proof. For convenience, we denote /(r) = i~^F^^{ir). Clearly, from the above remark, 
/: M — )■ [— oo, +oo] is real analytic, with the exceptions of a possible number of points 
which are poles of F^^. It is an easy exercise to observe that all poles of must be 
simple: indeed, otherwise it would follow from equation (13. 3p that is identically equal 
to zero. Equation (13. 6p is re-written for / as 

(3.7) f{r) = f{rf-rf{r)-c, r G R. 

In particular, F'^^{ir) = f'{r) G M. It is known [1] that Xmiy^j^oo F'^^{iy) = 1, and /(r) = 
9F^^(zr) > r for all r G (0, +oo). We claim that in fact this must hold for all r G M, fact 
which excludes the existence of poles on the imaginary axis for F^^. Indeed, continuity 
of / requires that for this inequality to be reversed, there must be a point |/ G M so that 
f{y) = y. Then, from (13.71) . it follows that f'{y) = — c G (0, 1). But for the real analytic 
/(r) to cross below the first bisector as r moves towards y from the right, we clearly must 
have that f'{y) > 1. Contradiction. This implies also that there is no real point at which 
the limit of / from the right is — oo. That the limit from the right cannot be +oo at any 
real point is trivial: that would make, according to (13.71) . /' tend to +oo, instead of — oo, 
at the same point when r approaches the point from the right. This proves half of (2) and 
all of (3). Next, the critical points of /: First, f'{s) = is equivalent to 

, s± a/s2 + 4c 

/(S) = n • 
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Since f{s) must be real, this excludes points s e (— 2-\/^, 2-^/^) (recall that c e (—1, 0)). 
Assuming such a point s > exists, we must have, from part (3), f{s) > s. However, 
that is clearly impossible, since c < 0. Thus, only negative s are possible, and for such an 
s < — 2\/^, we have 

s±Vs^ + 4c 

/(s) = ^ < o> 

or, differently stated, both any critical point and any critical value of / must be negative. 
We shall establish that indeed there exists a unique such s (depending of course on c), but 
in order to do that, we will first prove part (1) of the lemma. We claim that there exists 
— oo < go < so that / maps [go, +oo) onto [0, +oo). Indeed, otherwise /(M) C (0, +oo) 
and, as seen above, then /'(r) > for all r e M. Lagrange's theorem would require that 
there exists a sequence {r„}„ which tends to minus infinity so that /'(r„) — > 0. But this is 
impossible, since it would require 

lim /(rn)(/(r„) - r„) = c < 0, 

n— >oo 

while both terms of the product are nonnegative. So G /(M). Choose the largest point 
in f~^{{0}) to be go- Now it is clear that in addition / maps [go, +oo) bijectively onto 
[0, +oo); since /'(go) — —c> 0, it follows that, first, f~^{{0}) = {go}, and second, that the 
bijcctive correspondence extends to a strictly larger interval. We show next that this larger 
interval cannot be M. To do this, first let us assume towards contradiction that f'{r) > 
for all r e M. Then, of course, limr^._oo /('") = d exists and belongs to [— oo,0). We first 
assume that d e (— oo, 0). Then, by Lagrange's theorem again, we must be able to find a 
sequence {r„}„ which tends to minus infinity so that /'(r„) — > 0. As seen above, then 

lim f{r-n){f{r-n) - r„) = c < 0, 

n— >-oo 

or, differently said, either 

lim /(r„) = 0, or lim (/(r„) - r„) = 0. 

The first case cannot happen, since it would require that there exists a critical point of / 
in (— oo, go), which we assumed not to happen. The second case would require that 

lim /(r„) = -oo, 

n->oo 

which would contradict d G (— oo, 0) again. We consider then the situation in which 
d = — oo. To fulfill this condition, and in addition to avoid that f'{s) — for some 
s e (— oo, go), it is necessary that 



s r — Vr^ + 4c 

r < jyn < , 

for all r < —2^/^. (The necessity of the first inequality was proved before.) However, 
recall that /(go) = =^ /'(go) = — c e (0,1), so there must be then some point 
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t G (—00, go) so that f'{t) = 1. This imphes 

, . _ t ± y/t^ + 4(c + 1) r > if choosing + 
2 < t if choosing — 

which is an obvious contradiction. This, in addition, forbids the case of f'{t) = 1 for 
any t G (— oo,0]. Thus, there must be a critical point s G (— 00, — 2-\/^) of /. Since 
by equation (13. 7p f'{s) = =^ f"{s) = —f{s) > 0, any critical value of / is a local 
minimum, hence there exists only one such s, and f{s) is the global minimum of /. So 
/(r) G (/(s),0) for all r G (— oo,s). The previous arguments about the behaviour of / 
near —00 can be easily reapplied to show that lim,._^_oo f{r) = limr_5._oo f'i^) = 0. Finally, 
from part (4) above, it follows that f'{r) < 1 for all r G {—00, 0]. Thus, f'{t) > 1 for some 
t G M implies first that t > and second, that there exists a to > so that /'(to) = 1- 
Differentiating in (13. 7p and using part (3) gives 

/"(to) = /(to) - to > 0. 

Thus, /' increases at to whenever /'(to) = 1. In particular, this can happen for only one 
t G M, and thus /'(r) must tend to one with values strictly greater than one when r — )■ +00. 
But this contradicts part (3) (namely that /(r) > r for all r G M.) Thus it is impossible 
to have /'(to) = 1. This proves (5) and concludes our proof. □ 

The following lemma is trivial: 

Lemma 3.7. For a fixed c G (— l,+oo) there exists a t > depending on c so that 
F^,(C+) DR + it. 

Proof. Since F^^ is analytic, hence open, on the upper half-plane, and it increases the 
imaginary part, it is clearly enough to show that (i) lim^_^-|-oo ^F^d^ + i) = +oo and (ii) 
there exist M > N e [0, +00) so that < QF^^^x + i) < M for all x e R. To prove 

(i) , just observe that, since /ic is symmetric and has all moments, there exists a positive 
measure Ac, also having all moments, so that Xc(R) = J^u^dfic{u) and 

(3.8) F,dz) = z-Gxdz), ZGC+. 

Since Ac is a finite positive measure, hm^-^-i-oo Gx^x + i) = 0, so (i) follows trivially. Part 

(ii) is equally simple. We have that 

1 < QF^^x + ^) = 1 + / — — ^ r^dKiu) < 1 + Ac(M) = 1 + / u'^dfx.iu). 

Jm 1 + (a; - u)^ 

So the lemma is true for any t > 1 + v^d^^u). □ 

It follows now easily that F^^ is injective on the upper half-plane and F^^(C^) contains 
C+ + it for some t > depending on c G (—1, 0]. However, in order to prove our theorem, 
we need to find a larger set C D C+ which is mapped by F^^ bijectively onto the upper 
half-plane. To do this, we will show that for any t > there exists a D C"*" -|- it so that 
Ffj,^{Ct) = + it and F^^ is injective on Ct- This will clearly guarantee that has an 
analytic extension to + it for any t > 0, and hence (based on the previous lemma) a 
unique extension to C'^, concluding the proof of Theorem 13. II by an application of Theorem 
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13.21 Our strategy will be as follows: for a fixed t > there exists, by Lemma I3.6[ a unique 
s > go so that t = On the other hand, equation fl3.8p guarantees that there must 

be a number = N{t,c) > so that '^F^^{x + it/2) < t for all x G M, > N. Since, 
by Remark 13 ■4[ F^^ is locally injective around all these points, we conclude that there are 
three simple paths, one around s and two in {z G C''' : 'Rz > N} and {z G C"*" : 3^2; < — A^}, 
respectively, which are mapped by F^^ in M + it. We shall argue that these paths can be 
extended to a simple path pt containing all of them, with the property that F^^{pt) = M+^t. 
The correspondence, if existing, must be bijective, by Remark 13.41 and we will define Ct 
to be the simply connected component of C \ which contains numbers with arbitrarily 
large imaginary part. It will then be easy to prove that Ct has the desired properties for 
all t > 0. 

Lemma 3.8. With the above notations, there exists exactly one simple curve p = pt, 
symmetric with respect to the imaginary axis, passing through the point is G iM and so that 
^i^ciPt) = K + zt. Moreover, F^^ maps pt bijectively onto M + it. 

Proof. For an arbitrary t > 0, it is a consequence of Remark 13.41 and Lemma [3.61 that F^^ 
is conformal on a small enough ball centered at is. Thus, there exists a simple path p^, 
symmetric with respect to the imaginary axis, which is mapped bijectively by F^^ onto 
some interval {—e,e) + it for e > small enough. We shall show that p = pi extends 
analytically to an infinite path, denoted by pt, sent by in R + it, on which has 
no critical points and so that the limit at infinity of F^^ along either half of pt is infinite. 
Since is meromorphic, hence open, this will suffice to prove our lemma. Indeed, let us 
consider the connected component pf of F^^([0, +oo) + it) which contains is. It is clear 
that, as F^^ is meromorphic on C, the path pf must end either at infinity or at a pole 
of F^^, call it C- If it ends at a pole, it follows easily that F^^{pf) = [0, +oo) + it and 
the correspondence (by Remark [3T^ is bijectiveE Let us consider the second case, namely 
when pf ends at infinity. In this case, the possibility of having F^^{pf) = [0,d) + it for 
some < d < +oo must be discarded first: this would correspond to when F^^ has d + it 
as an asymptotic value at infinity along p^. Thus, let us show that 



Assume towards contradiction that this limit is finite, and call it x (the case when the 
limit does not exist is easily discarded). Of course, '^x = t > 0, 3ftx = d > 0. We shall use 
Equation (13.61) to obtain a contradiction: it follows from it that the differential equation 
satisfied by the inverse F~^^ (defined on a neighbourhood of [0, d) + it and with values in 
a neighbourhood of pf) is 



lim ^F^^{z) = oo. 



(3.9) 



1 



—w 



^ + wF ^^{w) — c 



This situation will turn out later not to occur, but at this moment this is not important for our proof. 
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Define r(t>) = ^F^^^ {v + it) , l{v) = QF^^^{v + it), < v < d. Equation (13.91) translates into 

(3 10) r'(v) = v{riv)-v)+t{t-iiv))-c 

[vir{v)~v)+t{t-L{v))-cf + [t{riv)-2v) + VL{v)f 

(3 11) l(v) = t(r(i;) - 2v) + vl{v) 

[v{r{v)-v) + t{t-t{v))-cf + [t{r{v)-2v)+vt{v)f 

As noted before, limy^a F~^^{v + it) = oo, so at least one of r(t>), l{v) must be unbounded. 
Thus, at least one of r'{v), i'{v) must be unbounded. From equations fl3.10p and (13. lip it 
follows that in order for any of r'{v), l'{v) to be unbounded, it is necessary that there exists 
a sequence {f „}neN C [0, d) tending to d so that 

lim [Vn{r{Vn) - Vn) + t{t - L{Vn)) - cf + [t(r(w„) - 2Vn) + Vni{Vn)f = 0- 



n—^oo 



As noted in the comments preceding Remark 13. 4[ L{vn) < t. Let us choose a subse- 
quence of {fnjneN, also denoted by f„, on which t{vn) converges. If it converges to 
i G (— oo,t], then we know that r(t>„) — )■ oo (since F^^(t>„) — )■ oo). But then, in or- 
der for the above displayed limit to hold, it is also necessary that L{vn) tend to infinity, a 
contradiction. So we must have that both r(f„) and i{vn) tend to infinity (plus or minus). 
Then lim„^oo K(r(f„) - Vn) + t{t - L{vn)) - cf = Q =^ lim^^oo ^nxj^ = t, or equiv- 
alently lim„_ ^ = i. Similarly, lim„_ = lim.^^ = -f. So f = -i, 

which implies d'^ = —t"^, a contradiction. Thus, it is impossible for F^^^{vn) to tend to 
infinity when f„ — )■ d, so F^^ cannot have a finite asymptotic value along pf. Since /ic is 
symmetric, this concludes the proof of our lemma. □ 

By proving the previous lemma, we have also proved that the inverse F~^^ of F^^ admits 
an analytic extension around i{0, +00) and around R + it for any t > 0. We shall argue 
that all these extensions agree with each other, and provide us with an analytic map F~^^ : 
C+ — > C which (1) decreases the imaginary part, and (2) satisfies limy_5.+oo F~^{iy)/iy = 1. 
Let us denote 

to = inf{t > 0|M + ir C F^,(C+ U M)Vr > t}, 
and So the unique number greater than go so that F^^(iso) = ito. As noted after the proof 
of Lemma [3.71 it is clear that F~^^ : C"'' + ito — )■ C"*" satisfies both (1) and (2). Clearly, this 
function has, by Remark 13.41 and Lemma 13. 6[ a unique analytic continuation to a small 
enough neighbourhood of ^[0, +00) in C^, which we will still denote by F~\ Now, for an 
arbitrary t G (0, to), we have proved in Lemma 13^ that F~^^ admits an analytic continuation 
to a small enough neighbourhood of R -|- which coincides on a neighbourhood of it with 
F~\ Since C"*" is simply connected and C"*" = (C"*" + ito) U \^ {R + it), we conclude 

te(o,to) 

that F~^^ admits a unique extension to the upper half-plane, with values in C. Let us 
denote C = -F^^(C^). It follows easily from the identity principle for analytic functions 
that F-\Ff,^{z)) = z for all 2 G C and F^^(F-i(z)) = z for all 2 G C+. Moreover, 
^F^XA > for all 2 G C+ U M, and if z G C", F^XA ^ C+, then it is obvious that 
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the inequality '^F^^{z) > '^z still holds. Thus, we conclude that both (1) and (2) are 
satisfied by F~\ Now our theorem follows: we define (p^^{z) = F^^{z) — z, z E C+. This 
function is obviously well defined, and maps the upper half-plane into C^, thus satisfying 
the requirements of Theorem 13.21 □ 

Since the classical Gaussian 7 equals /iq, we are now able to conclude its free infinite 
divisibility from Theorem 13.11 and Theorem 13.21 

Corollary 3.9. The classical normal distribution d'y(t) = -^e'^^^^dt is freely infinitely 
divisible. 

Next we shall discuss some properties of the Askey-Wimp-Kerov distributions fic with 
parameter c G (0,1). As mentioned in the introduction, numerical computations show 
that fic is not freely infinitely divisible for certain values of c > 0. More specific, direct 
(computer assisted) calculation of the Hankel determinants of the free cumulants of /ig/io, 
for example, shows that the 97*^^ determinant is negative, and thus the cumulant series 
of /ig/io is not the moment sequence of a positive measure on R; Theorem 13.21 allows us 
then to conclude that /xg/io is not freely infinitely divisible. (More such computations have 
been performed, and they seem to indicate that the size of the first Hankel matrix whose 
determinant is negative rather tends to decrease as c > increases: for example, the 83*^ 
Hankel determinant corresponding to Hi is negative.) 

However, the family {/Xc : c G (0,1]} turns out to be of some interest from the point 
of view of the arithmetic properties of free additive convolution. This subject is not new 
(implicit results on the arithmetic of free convolutions can be found in many works), but it 
is the rather recent preprint |20] that has first addressed the problem of the decomposability 
of measures in "free convolution factors" in an explicit and systematic way. However, the 
subject is by no means exhausted, the number of results is rather small (we would like 
to mention among them a remarkable idecomposability result given in p]), so we feel it 
is worth mentioning the following by-product of equation (13. ip and our main free infinite 
divisibility result from Theorem 13.11 

Indeed, it follows from the continued fraction expansion (13. ip and analytic continuation 
that for any c G (—1, 0], 

(3.12) G,^,,{z) = '~f^'f\ zee. 

In addition, for any fixed c G (—1,0], the dilation transformation fic p-l induced by 
Ffj^i{z) = ^f^ F^^{zy/c + 1) provides us with a probability fil of variance one, and thus 

(3.13) F^i{z) = z - v/^G^^^,(;2v/FTl) = z - G-^^^,{z), z e C, 

where of course flc+i is a probability measure obtained by a dilation with a factor of 
\/r+~c of fic+i- It is noted in [3 Theorems 1.2 and 1.6] (see also references therein) that a 
probability measure A with variance one and first moment zero is freely infinitely divisible 
if and only if there exist two probabilities u, p onM. so that 

(a) F^{z) = z- G^{z) for all ;z G C+, 
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(b) Fx{z) = z — Gussiz), z E C+, where A = (p^^) and S is the centered Wigner 
(semicircular) distribution of variance one. (Operation l±) is called Boolean convo- 
lution - see [39].) 

We apply this observation to A = /ij to conclude from fl3.13p that the Voiculescu transform 
0^1 of fil is also the Cauchy-Stieltjes transform of a probability measure Uc playing the role 
of i> in (a) and (b) above and moreover 

Pc+i = t'c ffl '5. 

This provides us with another interesting decomposition result in the arithmetic of free 
additive convolution, stating that 

Remark 3.10. For any c G (—1,0], the Askey-Wimp-Kerov distribution /ic+i can be 
written as a free additive convolution of the Wigner law with another probability Uc on M. 
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